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Abstract 

A review of light-cone, covariant and gauge-invariant wave-functions of mesons is pre- 
sented. They are basic non-perturbative objects needed for hard exclusive processes and 
for the method of light-cone QCD sum rules. The emphasis is on the vector mesons and 
a new (model-independent) way of computing the mass corrections of vector-meson wave 
functions is given. 
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1 Introduction 



One of the main goal in QCD is to understand the non-perturbative aspects of this theory. 
In this report, I describe a way of parameterizing these non-perturbative effects in a co- 
variant and gauge-invariant way : the "meson wave functions" or "hght-cone distributions 
amphtudes of meson" , with the emphasis on vector meson (the pseudoscalar case, hke the 
TT-meson, is simpler and the results can be found in 0). By definition, a "wave function" is 
the matrix element of (gluon and) quark operators on the light-cone between the vacuum 
and a light meson. It is conceptually different to the OPE applied to QCD (like in the 
method of QCD sum rules [Q) where the non-perturbative inputs are vacuum expectation 
values of local operators. 

The standard approach to wave function (due to Brodsky and Lepage |^) considers 
the parton decomposition in the infinite momentum frame. A mathematically equivalent 
formalism is the light-cone quantization]^. Here, the concept of wave function is defined in 
a different way. Although it is less intuitive, it keeps the Lorentz and gauge invariance, in 
order to use the equation of motion of the QCD. With this approach, one gets directly the 
non-perturbative objects needed for the method of light-cone QCD sum rules|2T|, which 
allows to compute exclusive decays of heavy hadron beyond perturbative QCD. 

To compute a whole non-perturbative function seems an impossible task. Therefore 
the "conformal expansion" is used (it is like a partial- wave expansion, but referred to the 
conformal invariance of QCD), and the unknown input come from local matrix elements, 
which can be computed with the help of QCD sum rules[Q. Within this expansion the 
renormalization-scale dependence is also under control. 

Like in DIS, the different wave functions are classified by their increasing twist. In 
order to be more precise when using these wave functions, the non-leading twist ones are 
needed. A part of these non-leading twist contributions comes from the mass correction, 
which can be important if K* is considered for example. 

In this work, after two sections of preliminaries and definitions, I present the conformal 
expansion of these wave functions : I introduce the conformal group, the notion of confor- 
mal operators on the light-cone, the link between conformal invariance and renormalization, 
the matrix elements of conformal operators. Then I give an example of computation of 
non-leading twist three-point wave function. Finally, I give methods of computing mass 
corrections, with some new formulas which gives the mass corrections directly from the 
leading twist wave function (without making a conformal expansion). I add some appen- 
dices which contain lists of wave functions and technical results. 

2 General framework 

The two basic objects are (for the p-meson) : 

1) The two-point wave functions, extracted from this kind of matrix element : 




(2.1) 



2) The three-point wave functions, extracted from this kind of matrix element : 




(2.2) 
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Where x is almost on the hght cone, w G [0, 1], F any kind of product of 7^ matrices and 
[x, y] a path-ordered gauge factor along the straight line connecting x and y : 



[x, y] = P exp 



1 

ig / dt{x-y)^A^'{tx + {l-t)y) 




(2.3) 



This gauge factor is a way to introduce interaction (see [|22|), it ensure the gauge 
invariance of these non-local matrix elements (I will omit to write it sometimes). 

These matrix elements depend on three vectors : 
- : momentum of the p-meson 



- 

X I, 



(A) 



polarization vector of the p-meson 



With the relations : 

(2.4) 
(2.5) 
(2.6) 

The parameterization of these different matrix elements is based on the Operator Prod- 
uct Expansion on the light-cone^]. So one needs light-like vectors p and z, built from P, 
e^'^^ and x : 





p2 


= 


e(A) 


e(A) 


= -1 


P 


e{A) 


= 



p 








such that p ^ P in the limit — > and z x for x"^ —>■ (see 



^ 2^pz- 



xP — \ {xP) — x'^m'i 



Useful scalar products are 



zP = zp 



xP) — x'^m^p 
2 



Ipz 



(2.7) 
(2.8) 



(2.9) 
(2.10) 

(2.11) 
(2.12) 



The polarization vector e*-^-* can be decomposed in projections onto the two light-like 
vectors and the orthogonal plane : 



,{A) . 



Z) 



Note that 



pz 

(eW ■ z 



2pz^'' 



Zu \ +e 



(A) 



(2.13) 
(2.14) 
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For a practical use, one defines the projector onto the direction orthogonal to p,z : 



and these notations : 



at = 9tlU {Pt^Zy + PuZf,), 



pz 



(2.15) 



a* = a^p^ / {pz) 



(2.16) 



3 Definitions of twist 



The different contributions to the matrix elements |2.1| , p.2| are classified by their "twist", 
but there are two different definitions of this concept : 

For local operators, the "twist" means dimension minus Lorentz-spin. To extract the 



contribution of a definite twist for non-local matrix elements like |2.1j , |2.2| , an expansion in 
local operator (OPE) has to be done and all local operators which have a definite twist in 
the sense given above must be re-summed. This will give the non-local term which has a 
definite power in x^. To understand how this definition of "twist" means "power in x^", 
see 0. From now on, this definition will be called "theoretical twist". 

The second definition is built in analogy to partons distributions. In the latter case, the 
different structures of non-local matrix elements are separated with their different powers of 
Q (the hard momentum transfer) in deep inelastic scattering : a term of twist t contributes 
to the inclusive cross section with coefficients which contain t — 2 or more power of 1/Q. 
A good description of this classification can be found in . For the matrix elements |2.1| , 
2^ , the same type of definition can be used, like in ^ : structures with p ■ z, with p^ or 
with e*^'^-' ■ z gives one power of Q, but e^f^ behaves as Q". From now on, this definition will 
be called "physical twist". 

As an example, one gives the parameterization of a two-point matrix element : 



{0\u{z)j,[z,-z]d{-z)\p-{P,X)) 



fp^p 



P-z Jo 



,(A) 



1 e(^) ■ z 

-Za- T^mJ / due'^P-'g3{u,iJ.^) 



2 '-{p-z] 







where 



(3.1) 



e = 2M-l, (3.2) 
fi is the renormalization scale and fp is defined by the following local matrix element : 

(0|n(0)7^rf(0)|p-(P,A)) = /pm.eW. (3.3) 



All the three functions (f)\\,g^^\g3 are normalized like this : 



du(j)\i{u) = 1 



(3.4) 
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Following the classification of 0|| is a ("physical"-)twist 2 contribution, g^^^ is 
(" physical" -)twist 3 and is (" physical" -)twist 4. But when one looks at the defini- 
tion given above, one would argue that these contributions should be twist 1, 2, 3. In fact, 
the authors of |p used another definition of twist, which just shifts it to one unity, based 
on the light-cone quantization formalism. The parameterization of all matrix elements 2^ 
and |2.1| up to twist four can be found in the first appendix. 



A physical meaning for the different functions (f)\\,g^^\g3 can be given : they describe 
the probability amplitude to find the p in a state with a quark and an antiquark which 
carry momentum fractions u for the quark and 1 — u for the antiquark respectively and a 
small transverse separation of order 

Roughly speaking, one uses the "physical twist" to define the different wave functions. 
But in order to compute them, an expansion in "theoretical twist" is made. 



4 Conformal expansion 

When the different wave functions which can be extracted from non-local matrix elements 



like ^]T] and ^]2| has been classified, one needs to compute them. The first idea is to 
expand them in the different moments, like in But it is like doing an expansion in local 
operators (OPE). This must be avoided if these wave functions are needed for a light-cone 



QCD sum rule (see [21 



The method used to compute wave functions is the " conformal expansion" : the idea is 
to expand the wave function in a series of polynomials, whose coefficients are renormalized 
(almost) multiplicatively. Here is an example for 0|| (with ^ = 2^ — 1) : 



</>„(n) = QuuJ24cT{0,^ = 2« - 1 (4.1) 

n=0 

where Cn'^ is the Gegenbauer polynomial (see appendix |CD and all are coefficients which 
are renormalized multiplicatively (to one loop). These ai can be computed using the usual 
QCD sum rules and exact relations (see section ||). 

In this section, the method of building such an expansion will be described. For it, 
one needs to know what is the structure of the conformal group, in what sense is QCD 
conformal invariant, how to build local and polylocal conformal operators and finally what 
happens for the matrix element of polylocal conformal operators. 

4.1 Structure of the conformal group 

The conformal group is a Lie group which contains the Poincare group and the dilatation. 
To describe its structure, it is simpler to work directly with the algebra. The generators of 
the conformal group are the generators of the Poincare algebra and M^,y, the dilatation 
D and of the special conformal transformations K^. 



They satisfy the following commutation relations]^, |T^, |T^ : 

[D,K^] = iK^ (4.2) 
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[K^,PJ\ = -2i{g^,D + M^,) (4.3) 

[D,M^,] = (4.4) 

[D,P,] = -iP, (4.5) 

[Kp,M^,] = i{gp^K,-gp,K^) (4.6) 

[K,,K,] = (4.7) 

plus those of the Poincare algebra. The action of these generators on a field $(x) with an 
arbitrary spin and canonical dimension / is the following : 

[P„Hx)] = -td,<!>{x) (4.8) 

[M^,,^x)] = [tix,d^-x,A)-^,..Mx) (4.9) 

[D,^{x) = -i{x^d^ + l)^{x) (4.10) 

[K^,<!>{x)] = -z(2x^x% - x^^^ + 2V - 2^a;,S^,)$(x) (4.11) 

where S^j, is the generator of spin. It acts on a fermion field ip and a gauge field Gap in 
the following way : 

^f^ui^ = \(y^ui^ (4.12) 

^fivGap = i{g^iaGup - guaGf,^) - {a^ (3) (4.13) 



4.2 Conformal invariance and QCD 

Anyone can convince himself that the unrenormalized QCD-action (with massless quark) 
is conformal invariant. But this symmetry is broken by two effects : 

1) Renormalization 

2) Quark-mass terms 

But a careful analysis of these effects shows that the conformal expansion of wave 
functions like in [4.1| is still relevant. 



4.2.1 Renormalization effect on the conformal invariance 



Suppose one has a conformal operator in non-interacting and massless QCD : it is an op- 
erator which does not mix with other ones under conformal transformation (such local and 
polylocal operators are constructed in the next subsection). Following [|15|, one can prove 
that such operators do not mix under renormalization at the leading log approximation (if 
simple derivatives are transformed into covariant ones). It means that conformal symmetry 
(and conformal expansion for wave function) is still a "good" symmetry of massless QCD. 

Consider an operator which renormalizes multiplicatively. Its matrix elements F satis- 
fies the Callan-Symanzik equation : 



r(p,/^;^) = 7r(^)r(p,/i; 



(4.14) 



where p = {pj^, . . . ,p^) is a set of momenta and 7r is the anomalous dimension of the oper- 
ator considered. If this matrix element comes from a massless theory, it is an homogeneous 
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function of p and jj,. The Euler's theorem yields 



■ d d 



(4.15) 



where dr is the canonical dimension of the operator considered. So equation [4.14| can be 
transformed in 



r(pp, /i; g) = [dr - 7r(^)] r(pp, /i; g). 



(4.16) 



If one considers the theory at one loop level, l3{g) can be neglected because it starts 
with g^] one gets : 

d 

This last equation shows that F satisfies the dilatational Ward identity, so if the in- 
teraction is switched off, one gets a matrix element which transforms in itself under a 
dilatation : it is the matrix element of a conformal operator. 

So if an operator is renormalized multiplicatively at one loop order, it is conformal if 
the interaction is switched off. Is the converse true ? 

Suppose one has a set of operators. A basis which diagonalizes the anomalous dimension 
matrix at one loop order is built. From the preceding arguments, one knows that this basis 
is made of conformal operators. If one directly builds a basis of conformal operators (like 
in the next section), almost the same basis of multiplicatively renormalized operator is 
obtained, except if the eigenvalues of the anomalous matrix (or in the representation of 
the conformal group) are degenerate. 

So conformal operators renormalize almost multiplicatively at one loop order, they can 
only mix with the ones which have the same conformal representation. The coefficient of 
the dilatational Ward identity transforms like this 



dr ^ dr- 5'%2 



(4.18) 



when the interaction is switched on at one loop order. (7^72 is the value of the anomalous 
dimension 7r(5') at the order g"^ : 



ir{g) = g'^i2 + ^74 + • • • 



(4.19) 



This property of conformal operators works also in the leading logarithmic approxima- 
tion, otherwise it would contradict this property at one loop order (see Wl 



In the anomalous dimension of some conformal operators is computed and for 
these operators, it is checked that they renormalized multiplicatively to one loop. 

Some complications arise when one deals with gauge theories. In order to have gauge- 
invariant operators, usual derivatives have to be replaced by covariant ones. But in that 
case, it is not sure that the diagonalization of the anomalous dimensions matrix is really 
possible. But one hopes that for color-singlet operators, it works (see |]12|)- The other prob- 
lem is that an operator which transforms in itself under a dilatation may not be a conformal 
one (the action of special conformal transformations iiT^ may introduce some gauge- variant 
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terms). But it seems very hard to construct such " almost" -conformal operators, so it is 



assumed that they do not exist (see [p 



The coefficients of the conformal expansion of a wave function (like the a„ in equation 
|4.1| ) are matrix elements of local conformal operators (see sections |4.3.1| and ^). So they 
are renormalized multiplicatively, that is why conformal expansion like in ^?T] make sense. 



4.2.2 Quark-mass effect on the conformal invariance 

Masses of light-quarks break explicitly the conformal invariance of the QCD-Lagrangian. 
So they can break the multiplicative renormalizability of conformal operators. But if the 
renormalization scale is high enough compare to the light quark masses, this effect is small 
and the conformal expansion of wave functions is still a "good" expansion. 

Quark masses can also appear explicitly in the twist-expansion of matrix elements of 
non-local operators. It that case, it introduce some new wave functions, but the latter can 
also be expanded in a conformal way. 

So the effect of quark masses can be controlled and the conformal expansion is still a 
good way to compute the wave functions. 



4.3 Construction of local and polylocal conformal operators 
4.3.1 Local operators on the light-cone 

For the use of QCD light-cone sum rules, one deals with fields ^{z) varying on the light 
cone (z^ = 0) or almost on the light-cone. But the representation of the conformal group 
is much simpler on the light cone, that is why a light-like basis z^, was introduced at 
the beginning of this chapter (equations |2.9| and |2.10|) . 



Following 10, one can see that only the components P., D, M*. and K^, of the algebra 
of the conformal group only act non-trivially on the field ^{z) (the meaning of . and 
J, is given at equation |2.16 ). These generators form the algebra of a subgroup called 



collinear conformal subgroup 50(2,1) = S'f/(1, 1) = SL2{R) which generates projective 
transformations on the line (see e. g. ||12|| ). 

This algebra can be brought into a more standard form, with the following linear 
combinations : 



(4.20) 
(4.21) 
(4.22) 
(4.23) 

JI-J3-J+J-. (4.24) 



J+ = Ji + iJ2 = 






J_ = Jl — iJ2 = 






J3 = 






E = 


>- 




J' = 


72 _ 72 


- 7^ 
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Then the commutation relations are 



[J+,J-] 
[J3, J±] 
[E, J,] 



-2J3 
±J± 






(4.25) 
(4.26) 
(4.27) 
(4.28) 



In order to build a representation of this subgroup, all the generators P., D, M^. and 
must reduce to differential operators. Equations |4.8| to [4.11| imply that fields must be 

eigenvectors of the operator S^,., that is fields having fixed projections (s) of the Lorentz 

spin on to the line : 

S,.$(z) = is^z) (4.29) 

A spinor field ifj has two components with the projection s = ±| : j.ip and 'j^ip- O^ig 
can see it using equation ^.12| : 



1 i 

^,.7.^^ = -cT,.7.?/' = -27-V' (4.30) 
1 i 

T.^.'j^i! = -a^.'j^ip = +-7*^ (4-31) 

For the gauge field G^^, there are three possible value of the projection : s = —1, 0, 1. 
The different components are (see equation [4.13|) : 

= +iG.^ (4.32) 
= -iG,± (4.33) 
S^.G,. = J:,.G±± = (4.34) 

where ± is a component in the plane orthogonal to the plane of the light-like vectors Zfj_ 
and Pfj^ (the projection operator associated to this plane is g^^,, defined by the equation 



If one has a field ^{uz) {u is a real number) which has a fixed projection s of the 
Lorentz spin on to the line 2^, the generators of the coUinear conformal subgroup act in 
the following way : 

\J+,^{nz)] = ^l^-fM (4.35) 
[J3,Huz)] = 1 + s + 2n^^ $(nz) (4.36) 
[J^,^{uz)] = V2(^u{l + s)+u^-^<^{uz) (4.37) 

[E,<!>{uz)] = ^(/-s)$(mz). (4.38) 

From this relation the irreducible representations of the collinear conformal subgroup 
5*0(2, 1) which contains the non-trivial conformal transformation on operators on the light- 
cone can be constructed. These representations are classified by the eigenvalues of the 



8 



Casimir operator J^. Equations [4.301 to [4.34| describe the construction of operators $(m2;) 
which are eigenvectors of : 



[j2,$M]=j(j-l)$M (4.39) 

where 

3 = \{l + s) (4.40) 

j is called "conformal spin". For $(uz), the algebra is reduced to differential operators, so 
one can work with functions of one real variable instead of quantum fields. Having in mind 
the application to the computation of wave functions, one can make the Fourier Transform 
of the quantum field $ : 

^{uz) = j c/ae-'"°P"0(a). (4.41) 

In this formalism, the irreducible representations of the collinear conformal subgroup 
can be built, classified by the eigenvalues of and — J3. They are the set of the following 
functions : 

1 / a 



Equations [4.35| to |4.38| give the action of the generators of the collinear conformal 



subgroup on the states \j, n) : 



J- li, n) 



{j + n)\j,n + 1) 
{n-j)\j,n-l) 
-n\j,n). 



(4.43) 
(4.44) 
(4.45) 



One can see that J+ transforms a state one step upper and J_ one step lower. The 
lowest value of n is j. 

To build a local conformal operator on the light-cone, one has to have fixed projection 
of the Lorentz spin on to the light-cone (equations [4.3CI| to |4.34| ). If the operator is not taken 
at 2; = 0, one has to expand it in the states [4.42| . But in order to compute wave-functions. 



the 2;-dependence of polylocal operators must be controlled, like those in |2]^ and pT2. 



4.3.2 Polylocal operators on the light-cone 

Suppose that $i(mi2), . . . , $„(m„2) are n local operators which have a fixed conformal 
spin. What happens for the product ^i{uiz) . . . $„(m„2;) ? When one thinks of irreducible 
representation of the collinear conformal subgroup, the corresponding " Clebsch-Gordon 
coefficients" of a tensor product of different irreducible representations are needed : 

\j,n) = ^rni,...,..,njji,ni) . . . |j.,n,) (4.46) 

ni + ...+7ife=7i 

This seems a very complicated task, but one has to recall what is the conformal collinear 
subgroup : it is the Lorentz group in 2 -|- 1 dimensions (50(2,1)). So the irreducible 
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representations | j, n) can be interpreted as relativistic particles in an abstract or internal 
space of 2 + 1 dimensions. With this point of view, the state | j, n) is a particle with a mass 
j and an energy n (that is why n > j). 

Consider a system of k particles. The lowest invariant mass of this system is the sum 
of the masses of all particles, j = ji + ■ ■ ■ + jk- If the state of the lowest invariant mass 
has also the lowest energy, all the particles must be a their lowest energy level, rii = jj. 
This state is non-degenerate. In the language of irreducible representations of S0{2, 1), it 
gives : 



k 



^ ji, nmin = jmin ) ~ \jl,jl) ■ ■ ■ \jk,jk) (4.47) 

i=l I 

al''-^ . . . al'"-' (4.48) 



If n > j, the operator J+ = X]i=i can be used to "raise" the "energy" of the state. 
Since J+0(a) = -^^a(j){a), one gets : 




jmin = y^3i,n) ~ («i + ... + «fc)" (4.49) 
~ (ai + . . . + afc)"- . . «f (4.50) 



To have the whole sum for the state \j,n) (equation [4.46|) , not only the states at rests 
are needed but also the higher ones {\ji,ni) with rii > ji). In that case, there are 
different possibilities and the Clebsch-Gordan coefficients contains binomial coefficients 
(see |]12|). For a bilocal operator, one gets : 



(4.51) 



\j,n) = (ai + as)" ^ ^ \ ni^^ ) l^^'-^^ + ^i > \j2,j2 + n2 > 



ni+n2=j-ji-j2 



The summation gives : 

\j, n >~ {a, + a2)"+^'+'(l + - 0'^'^"'^^T4f "'^0 (4-52) 

with 

i = (4.53) 

where Pn^''''^\Cj are the Jacobi polynomials [|T^] (see appendix 0). Analogously, for product 
of three local operators, one gets Appell polynomials ||16[j (see appendix |Cp. The important 
thing to remark about equation [4.52| is that the only dependence in n is in the factor 
(«! + a2Y^^^^, the other parts of this equation depend only of the conformal spins ji, j2 
and j. 

An important property will be needed for the next subsection. For any kind of irre- 
ducible representation of the collinear conformal subgroup | j, n) which is constructed from 
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the tensor product of fc- irreducible representations, there is the following homogeneity 
property : 

\j, n) = „(ai, . . . , afc) = (ai + . . . + afc)""^-'"Vi,n(ai, . . . , d^) (4.54) 

where 

rv ■ 

a, = (4.55) 

ai + . . . + ttk 

This property comes almost directly from the constructions of the Clebsch-Gordon co- 
efficients (equations |4.50| and |4.51| ). Hence irreducible representations \j,n) of the collinear 



conformal subgroup on the space of k variables induce the irreducible representations on 
the functions defined on the simplex ai + . . . + = 1. This property will be used for the 
next subsection. 



4.4 Matrix element of polylocal conformal operators 

Suppose that one has the matrix element of a product of local conformal operators 
$i(mi2;) . . . $fc(wfc) on the light-cone between the vacuum and a massless meson state h 
of momentum p^. This matrix element can be parameterized like this : 

(0|$i(Miz)...$fc(Mfcz)|%)) = j P«,e-^^^('^i°i+-+'^'="'=V(«i,...,«fc) (4.56) 

where 

Vai = I dai . . . dakS \ y — 1 \ (4-57) 




According to the property of equation [4.54|, the wave function 0(ai, . . . , a^) can be 



expanded in different parts which have a fixed conformal spin j. The minimum value of the 
conformal spin is the sum of the conformal spins of every conformal operators ^i{uiz). The 
different parts of this conformal expansion are mutually orthogonal polynomials (Jacobi 
polynomials if there are only two operators) and form the complete set of functions on the 
simplex ai + . . . -|- = 1. 

In front of each of these polynomials there is a coefficient which is renormalized multi- 
plicatively. This conformal expansion is justified by the construction of irreducible repre- 
sentations of the collinear conformal subgroup. 

The first term of such an expansion is called "asymptotic" wave function and can be 



easily computed (equation |4.50|) : 



0,3(«i,...,«,)- r(2j0...r(2j,) "1 •••"^ ^^-^^^ 

with 

^I)a0as(ai) = 1 (4.59) 



Now one can understand better the meaning of the example given at the beginning of 
this section : 



6uuJ24cT{0- (4.60) 



n=0 
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(f)\\{u) is the contribution of the matrix element {0\u{z)'y^[z, —z]d{—z)\p~ {P, X)) where 
each quark field has a positive spin projection s = +| (see equation [4.30| ). Knowing 
that a quark field has the canonical dimension / = |, each quark field has a conformal 
spin jg = 1 (see equation [4.4U| ). So the asymptotic distribution amplitude [4.58| equals 
0as(ctq, «g) = Qaga-g and has the conformal spin j = 2. The higher terms in 4.60 corresponds 
of higher values of j and n in equation [4.52| . Denoting u = ag with ag + a-g = 1, one gets 

^^^"^ are proportional to the Jacobi ones 
1, see ^321). 



the expansion [4.60| (the Gegenbauer polynomials C„ 
Pn'^^ (see Appendix C) which appear when ji = j2 



5 Example of computation 

In this section, the method of computing wave functions (defined in the appendix A) is 
illustrated with non-leading twist three-point distributions for the p-meson (following |]TI|). 
One starts with the definition of the chiral-even three-points distributions : 



{0\u{z)i'ja[z,vz]gGf,^{vz)[vz, -z]d{-z)\p (P, A)) 



fprrip 



p ■ z 

e(") ■ z 



^P{p-zf 
+ higher twists 



[Pt.9av-Pu9c,f,]^{v,pz) 

[Pf,Zu - PuZf,]pa'^{v,pz) 



(5.1) 



{0\u{z)-fal5[z,vz]gG^i.{vz)[vz, -z]d{-z)\p (P, A)) 
= ifpmp Pa[p^e';^l-p^e';^l]A{v,pz) 



p ■ z 
■ z 



^ {p ■ z)'^ 
+ higher twists 



[Pf^9L-Pugap]Hv,pz) 

[PfiZu - PuZf,]pa^{v,pz) 



(5.2) 



One wants to compute the non-leading twist distributions ^{v,pz), \E'(f,pz), ^{v,pz) 
and \E'(f,pz). The leading twist distributions have the (following) conformal expansion 
(see H) : 



A{v,pz) = J Pae'P^(""-"''+''"«)A(a) 



A{a) = SeOCsadanttq 



l+wt-i7ag-3) 



V{v,pz) = J r'ae^P^(""-"''+^"«)\/(a) 



(5.3) 
(5.4) 
(5.5) 
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(5.6) 



To compute $(f and the wave function of conformal operators is needed 



{Q\u{z)gGf,y'^.'^c,lbl*d{-z)\p{P,\)) = fpinl 



{Q\u{z)gG^u-i^-ial^-i.d{-z)\p{P,\)) = fpml 



,(A) 



z 



p ■ z 

e(^) ■ z 
p ■ z 



[PlJ^aiu - Pu9at^] (5-7) 
[p^^9au-Pugt]'^^^('^^P^) (5-8) 



(this parameterization corresponds to the highest conformal spin for Gf^,^, although this 
gluon field does not look explicit like a conformal operator, see equations [4.32| to [4.34| ). 

These two matrix elements have fixed conformal spin, so a conformal expansion (with 
the help of the Appell polynomials) can be done : 



1 + ujl^oiag - 3/2a„) + ujo^ii'^g - 3"^) + • • • 



Ti 



= K^^QOaaal [l + cu|J,(a^ - 3a„) + uj^^^{ag - 3/2^^) + 
Then it is easy to get $ and $ from ^'^^ and : 

$(«) = -\^^Ha) + ^^Ua) 



1 r$n, 

2 L 



$(a) = - [$^Ha) - $^^(a 



(5.9) 
(5.10) 

(5.11) 
(5.12) 



The G-parity (charge conjugation x isospin) applied to the matrix elements and 
implies the following relations : 



^1,0 — ^0,1 



"1,0 — ^0,1 

Finally, the expression for $ and $ (knowing that l — ag — au — ctd = 0) is the following : 
^a) = K^^30al{il-ag) + ^jj, [«,(1 - «,) - 3/2(a^ - a^)] 



Ti 



(5.13) 
(5.14) 

(5.15) 



+ ul\[ag{l- ag) -6auad] 



$(a) = K^^30al{au -ad) 1 + l/2iu[^Qi5ag - 3) + cul^^a 



,n 



(5.16) 



(5.17) 



Now, the different constants K^^, cjj q and ujI\ has to be extracted from matrix elements 
of local operators. 

To find K^^, the definition of ^{v,pz) in |5.2| is taken. It implies : 



fpm%e^^^ ■ z)g^, I Va<P{a) = {Q\uG.^^,^,d\p{P,\)) 



(5.18) 
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Using the equation 



Bl8| which defines (3 and (4, 



{0\ugG^^,-fa-f5d\p (P, A)) = fpiTipCs 



1 

3' 



(5.19) 



one obtains for K^^ 



-^C3 + ^C4 



(5.20) 



T I T I 

The two other constants uj{ q and uJq\ have to be computed. The idea is to relate these 
coefficients to local matrix elements. To do it, the matrix elements ^ are expanded to 
the order (zp). It gives 



(5.21) 



{0\ug D .gG.u + gC^ D . j 7-7a757*f^|0) = 
= ifpmlie^^^ ■z){p- z)gi^ j Va{au - ad)^^^{a) 

j Va{a^ - a,)$Ti(«) = K^i Q + l^n _L _ 



(5.22) 



{Q\ug [D .gG.,j 7-7a757*c?|0) 



fpm%e^''>-z){p-z)gt I Vaa,<l>^'{a 



(5.23) 
(5.24) 



These equations can be also rewritten like this (replacing K^^ by its value in |5.20| ) : 



(0|n -i D .gG.y + gG.^i D . i-iad\p{P, A)) 



-fpmy^lP.g, 



- I --Cs + -C4 I + —Ri - —R2 
6 V 3^ 3^ / 14 14 



+ . . . 



(5.25) 



u(i D .gG.A -fal5d\p{P,\)) 



fpmy^\P.g^, 



+ 



(5.26) 



where 



R2 ^ K^^cul^ 



(5.27) 
(5.28) 
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Because of the G-parity symmetry, one can suspect that only two local operators are 
involved (apart from operators with total derivative). One defines : 



(5.29) 
(5.30) 



Olpt^u = u(-D0 gG^^ + gG^^ B/j) '-fad 
The different Lorentz-structures for O^^^^ are the following : 

(0|O^^^Jp(P,A)) = 

-- "5 1 



= < e 



,(A) 



5 1 

—ml {Pag^v + P^gav) - 

5 1 

PaPpPu - {Pagpn + Ppgan) " ■^rnlP^.gap 



1 2 

m„ 

24 " 



+ ei-') (P^(7;3. - Pugp^.) D, + ej,"^ (P^^«, - P,g^^) E, 
- gap {pA'^ - PA'^) 



(5.31) 



And the constants A2 to F2 can be similarly defined, replacing by in this last 
equation. 

Traces can be taken in order to reduce the number of independent constants (it works 
for both and O^) : 

gap: B + C-D-E + 4:F^0 



ga^ : -e^u^Pp{B + 3C + P) - P,{D + 3E - F) = {0\O^^^MP, A)) 
gp^ : -eWp,(3P + C + F) - ei^^Pa{3D + E - F) ^ {0\Oa^^MP, A)) 
One defines a^, a^, and P in the following way : 

{^lOl'p^MP^ A)) = {ef^P, + e(")p^)ai,2/pm3 + higher twist 



(5.32) 

(5.33) 
(5.34) 



(5.36) 



(0|Ci'|^Jp(i", A)) = (eWp, + eWp«)6i,2/pm; + higher twist 

These higher (theoretical-)twist terms are not taken into account, because they are 
bigger than twist 4. This approximation gives : 



(5.37) 

(5.38) 
(5.39) 



B = 


D = 


1. 


- 36)mJ/, 


E = 
F = 


C = 



8^ 





The goal is to express the unknown coefficients a;o\ and a;| ^ in function of ^4 • • • P and 



a, b. 
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For it, the definition of (3 and ^4 in ID (equation p.8|) is used. It can be written like 
tliis : 

■ 1 . 1 



(0|m5(G.^7q75c/|p(P, A)) = e^^^g^^ufpfn 
Equation ^.311 implies : 



-3(3+30+... 



Ol:.i\p{P,X))=e'P.g, 



r/prnMi,2 + Ci,2 + ^1,2 + 



(5.40) 



(5.41) 



Taking |5.25| and |5.26| with the total derivative of |5.40| (with id. —>■ P.), combining [5.38 
with |5.41|, one gets : 



7 -T7C3 + 77C4 + TT^l - 77^2 = --[A 



14 



-77C3 + 7tC4 + —Ri + —R 



14 



14 



14 
1 



3ao 



C3 + 77C4 + -[Ai-bi + 3ai] 



(5.42) 



(5.43) 



The constants Ai^2 and Bi^2 can be extracted from the coefficients of the twist 3 wave 
functions A{a) and V{a). For that, the equations |5.1j and |5.2| are expanded at the order 



u 



D. i.gG., + gG.,^. D. ) d\p{P, A)) = f^m/^lp'-C^uJ^ 

7 ' 28 



(0|w7.75 D. gG.Jd\piP,X)) = -f,m/llp% ( ^ + -- 



The definition of A12 applied to these equations gives : 

'3 3 



Ai 
A2 



- + —^i 

7 28 ^ 



V 



(5.44) 
(5.45) 

(5.46) 
(5.47) 



Now, the three coefficients {K^^, Ri and R2 for $(a) and $(a) are expressed in functions 
of constants coming for the twist 3 wave functions and some unknown coefficients ai^2, ^1,2- 
Their number can be reduced using exact operator relations (with the equation of motion) : 

The equation 





implies (after taking the matrix element) 

hi = 
The equation 

implies 



0C3 — 7C4 
3 6 



iu-^alb [ G^a D( + G^a ) d 
^2 - ^2 = 0-1 - bi 



(5.48) 
(5.49) 

(5.50) 
(5.51) 
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And the equation 

Qdu0^al3 = 

= -12m7p ( Gpi3 Da - Da Gpp + {a ^)] d 



+ -d^ {daU-ffsd + dpU-iad) - 4:dpU [y/sGap + laGpp ) 75li 

8 ~ 8 ~ 

- ^dpu-iaGaalbd - -daUJaGap'j5d 

+ YSaisdpU'jaGapd 



implies 



24a2 



^ ^ ■((02)) + ^T^((Oo)) + ?C3-^C4 



^fprrip 



"ifpfnp 



with the following leading twist operator (symmetric and traceless) 

0^,af) = u (jpDaDp)^^^ d - traces 



(5.52) 
(5.53) 

(5.54) 



and {{On}) the different moments of the leading twist part of the two-point wave func- 



tion 



with 



{{On)) = Mr 



M^" = J duChi^^) 

{{{On}) is also defined by equation |b.l9[ ). 
Taking everything together, one gets : 

-i?2 = -(Ai + A2) + ai- — ^((02)) + ^ 



(5.55) 
(5.56) 



^ ((Oo)) + ^C4 



-Ri 



^ ^ ((Oo))-^C3 



36 



(5.57) 



(5.5J 



So the unknown coefficients i?i,2 (equation |5.27| and p.28| ) which are in the wave func- 
tions $ and $ are expressed in function of known coefficients, expect for ai which has to 
be computed by the method of QCD sum rules. 

With the same method, \1/ and \^ can be computed : 



^(a) 



A^* 1260 {adau)adauag 



^o + ^i(3a,-l) 



(5.59) 
(5.60) 



with 
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(5.61) 

^0 = ^C3 + ^C4 (5.62) 

So, in order to have the non- leading twist wave functions $, $, \1' and the coefficient 
tti has to be computed. But it is defined from the matrix element of a local operator in 
equation |5.35| , hence the method of two-point QCD sum rule described in |^ can be used. 



In |§] , there is a good presentation of the renormalization dependence of the coefficients 
appearing in the wave functions. The basic idea is to defined them from local conformal 
operators, where the renormalization dependence is better known (see section | 



6 Meson mass correction 



In this section, some new techniques (partially reproduced in |TT|]) to compute mass cor- 
rections of wave functions are described. 

The mass of the p-meson induce some non-leading (physical-) twist structure for the 
non-local matrix element. 

For example, for the chiral-even two-point distributions 



\u{x)-i^[x,-x\d{-x)\p (P, A)): 

e^^) ■ X 



+ 



P. 



>(A) 



P-X 

e^'^) • X 



1 1 



P. 



P -X 



due'^^-'' P(") 



2 2 

^p(^)| 



u 



1 

x, 

2 ^ 



X 



(P-x) 



u 



(6.1) 



some wave functions would disappear if the mass of the p were 0. Hence mass corrections 
are higher (physical-) twist contributions to non-local matrix elements (in the previous 
example, A^^^u) and C'^^\u) are twist 4, B^^^ is twist 5). In fact, the contribution 

can first be seen in the difference between the light-like basis (p^, and e^fl defined in 



2.101 , |2.9| and p.l3| ) needed for the conformal expansion and the vectors which characterize 
the p-meson (P^ and e^^''). These different basis introduce terms which are sort of 

"explicit" and "direct" corrections; but there are also the other corrections like A^^\u), 
P(^) and C(^)(m). 



6.1 Example of mass correction for a scalar theory 

Suppose that 0(x) is a scalar field with M a scalar particle of momentum P^ and mass m. 
If x^ is almost on the light-cone, we have the following expansion : 



(O|0(x)0(-x)|M(P)) 



g{u) + -m'^x'^g{u) + ... 



(6.2) 
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with 



(6.3) 



It is quite hard to compute directly g{u). In fact, the following approximation is done : 
in the Operator Product Expansion of 0(a;)0(— x), only the leading-twist local operators (it 
is the "theoretical twist", see section H) is taken. The leading (theoretical-)twist expansion 
of |S.2| is defined in the following way : 



(0| [ct>{x)<i^{~x\ , |M(P)) 



(6.4) 



In that case, f{u) can be computed in term of f{u). The first way to do it is the "brute 
force" method, i. e. doing a local expansion of p]4| : 



- Dx) (j) 



1. t. 



|m) = {iPxYMn + -m^x^n{n - l)M„_i(iPx)"-2 



(6.5) 



where 



dui^fiu) 



(6.6) 



The leading twist operator is always symmetrical in all his indices and traceless. 
It implies : 



J 1. t. 

2 " 

Pui ■ ■ ■ Pun ~^ T) ^ ^ 9uiUjPui ■ ■ ■ Pui-lPui + l ■ ■ ■ Puj-^Ppj + x ■ ■ ■ Pu„ + 0( 



m 



i<j 



(6.7) 



Hence 



I D X 
(Ps)" - 



1. t. 



|m) 



^^■(Px)"-2 + 0(m^) W(0„)) 



2n 2 

The comparison of this last equation with gives 

M„ 



nM„-2 



{{On)) 

((0„)) = M„ 



The integrated version of these equations is 



f{u) = 2{2u - l)9{u - 1/2) / dtf{t) + 2(1 - 2n)^(l/2 - u) / dtf{t) 



{6.i 



(6.9) 
(6.10) 



(6.11) 
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The same answer can be obtained using the techniques developed in [|T9l directly in an 
integrated form : 



(0| , \M{P)) = / du [e^^''% ^ /( 

^0 



u 



where 



m 



(6.12) 



(6.13) 



It can be checked that these last equations reproduce |6.11 

•1 nl 

) Jo 



2/" rfte*^^<2t~i) /■ dufiu)^^ [ dvv6{2t-l-^v) (6.14) 
Jo Jo Jo 



and 



^0 



dvvS{2t -1-^v) 



1 2/ — 1 

-^i-^ [e{t - i/2)e{u -t) + e{-t + i/2)eit - u)] 



(6.15) 



— (2t - 1) [9{t - l/2)9{u - t) - ^(1/2 - t)9{t - u)] 



(6.16) 

The goal is to reproduce this formalism in the physical case : a vector p-meson in QCD. 



6.2 Mass correction for two-point wave function : theoretical 
twist method 

6.2.1 Leading theoretical tvi^ist 

The idea is to compute the mass corrections in the leading theoretical twist approximation. 
In this part, the case with the two-point chiral even distributions is computed. 
The different wave functions are : 

{0\u{x)'y^[x, -x]d{-x)\p~ {P, \)) = 



fp^p 



P ■ X Jn \ 4 





, t\^ e(^) ■ X 



u, 



1 pW . r 



2^{P-x] 







(6.17) 
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The "brute force" method, i. e. the expansion of the distributions in moments (x 0) 
can be used. It gives : 



{0\uj,{zDxrd\p~-{P,X)) 



P. 



g(A) . 3, 
P-X 



1 1 

P^l^ '\\l\/rA 2 



^ ^ P-x 



(Px)"M, 



2 2 

a; m„ . „ o 



X 



~x I 



[P-x] 



^m%PxTM: 



(6.18) 



where the moments M^, M^, and are defined hke in 16. 61. 

The leading (theoretical-)twist of the operators in the left hand side can be written. 
The condition of symmetry and zero trace gives : 



{0\uj,{zD--xrd\p-{P,\)) 



g(A) . J, 



1 



p^i—-^iPxr + 

P ■ X n + 1 



e(A) _ p fill. 
^ ^ P- X 



+ 1 



.(A) 



P-X 



4 



^ IT, \n-2 - 1)^ 

- Px " 

(n + 1)"' 



1 



g(A) . 2; 



/n 2(^ - l)'^ 

-Xhtt; —{Px)ni'^ 



P ■ X)^ 



The comparison of these two equations implies : 



((On)) 
1 



n + l 
n — 1 



n-2 



{(On)) 
{{On)) 



m: 



c 



(n + l)2 

^"-'^^(o.)) 



((On)) 

(6.19) 

(6.20) 
(6.21) 

(6.22) 
(6.23) 
(6.24) 



(n + l)2 

There is another method to get the leading (theoretical-)twist part of the distributions 



6.17|, giving directly integrated equations between the different wave functions. The differ- 



ent relations for non-local operators in appendix are used : they give integrated forms 
for the leading twist part of two-point operators. 
One takes equation p.4| : 



\^{x)-fa^{-x)] 



sym. 







^ d - 

du— — il)iux)xil)i—ux) 

OXn 



(6.25) 



where x = x^'-j^. The subscript "sym." means that in the Operator Product Expansion 
(local expansion around x = 0), only the operators which have symmetrical Lorentz-indices 
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are taken. But for the leading theoretical twist, the traces has also to be removed, so the 
leading twist part of the right hand side of |6.25| is also taken : 



1 1. t. 



^ d 
du— — \%Ij(ux)xiIj(—ux)], 

OXr 



(6.26) 



D3D : 



The condition for the leading twist part of il){ux)xil){—ux) is the following (equation 

d 



dxn,dx° 



x)xip{—x)~\ J 







(6.27) 



Now the matrix element of ip{ux)xip{—ux) can be taken. The parameterization of the 
leading twist part up to corrections is the following : 



{Q\[u{x)x[x, —x\d{—x)\ J |p (-P, A)) 



(e^^) ■ x) / dug{u) 



2 2 
X Tfi 

e^^^^ + — ^/c(Pa;) + 0(mJ) 



where g{u) is the leading twist wave functions in the limit mp 0. 
The condition |6.27| can be applied; it gives (after some calculation) : 



(6.28) 



dug{u) [?>k{Px) + {Px)f^{Px) - ee'^""^] = 

Jo 



(6.29) 



The following ansatz for f^{Px) can be done (having in mind the case of scalar field in 
the previous subsection Ol) : 



Then equation |6.29| gives 



= a(0 / biv)e 



hiy) = 

This can be written in this form : 

{Q\[u{x)x[x,-x\d{-x\ Jp-(P,A)) 
dug{u) [(eW-x)e^«^^], ^ 



(6.30) 



(6.31) 
(6.32) 



(6.33) 



with 



[0 



Ji. t. 



(e(^) ■ x) 



4 Jo 



(6.34) 
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The equation |6.26| can be applied. After some algebra, one gets the distributions 
B'^^^u), C(^)(m), M^\u) and B^^\u) defined in^J^ in the leading (theoretical-) 
twist approximation : 



U 



9W 



[\uB'^'\u)e'^''A = [' dug{u) T rfte'*«^" 

Jo J 1. t. "'0 Jo 



(6.35) 
(6.36) 



1. t. 



1. t. 



dugiu) 



Jo 



dug{u)^'^ / dti 



duM'\u)e'^''A = f dug{u) [' dt^^ ['v^dve'"'^'''' 
J 1. t. Jo Jo Jo 



(6.37) 
(6.38) 

(6.39) 



There is no index "1. t." for A^^^ (u) because this distribution contains only the leading 
twist part of the matrix element (0|'u(x)7^[x, —x]d{—x)\p~{P, A)). So these last equations 
show that all the distributions can be computed when one has A^^\u), in the leading 
(theoretical-)twist approximation. They reproduced exactly the relations between mo- 
ments (equations |6.20| to |6.24| ). The advantage is that they are directly in an integrated 
form. 



6.2.2 Adding twist 3 and twist 4 terms 



The idea is to add contributions of (theoretical-)twist 3 and 4 without the gluonic contri- 
bution. In the appendix 0, there is the following operator relation ( p.lO|) : 



1 

nit. 4 



ip{x)xip{—x) = \^^J(x]x 

+ operators with gluons + . . . 
The matrix element of the term with a total derivative, gives 

1 



ip{vx)x%lj{—vx) 



(6.40) 



-X / dvv {Q\d u{vx)x[vx,—vx]d{—vx)\p{P,\)) 



m -X I dvv {Q\u{vx)x[vx^ —vx]d{—vx)\p{P,\)) 



(6.41) 



The 0(m°) term (equation |6.28|) in the parameterization of the matrix element in the 



right hand side implies : 



— -x^ / dvv'^ {{i\d'^u{vx)x[vx, —vx]d{—vx)\p{P,\)) 



2 2 /.I 

■ x) ^ 



(6.42) 
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So the leading twist part with the total derivative term is : 



(0| [u{x)x[x, -x]d{-x\ ^ 

dug{u) [(eW-x)e^«^^], ^ 



|P^(^,A)) 



tot. der. 



(6.43) 



with 



(e(^) ■ x) 



tot. der 



(6.44) 



This last form can be used to give some twist 4 corrections to equation |6.35| to |6.39| . 
Now the twist 3 part without gluon can be added. The matrix element of the equations 



D.7| is taken, suppressing the gluonic terms : 

(0| [M(a;)7^[x, -x\d{-x)\^ \p'{.P, A)) = 
= -^e^-a/j I dttx''{-iP''){0\u{tx)^p^^[tx,-tx\d{-tx)\p-{P,\)) 



(6.45) 



(0| [m(x)7^75[x, -x\d{-x)\^ \p (P, A)) 



-le^^o^p / dttx''{-iP''){Q\u{tx)^p[tx,-tx\d{-tx)\p-{P,X)) 



(6.46) 



Combining these two equations, replacing the right hand side of the second one with 
the parameterization of |6.17| , one gets : 



(0| [m(x)7^[x, -x\d{-x)\^ \p (P, A)) 



1 /•! 



{e^^\P ■ xf + P^{e^^^ ■ x){P ■ x) / du / dtt / dvv^B'^^\u)e''''^^'' 



^0 



+ ml f du f dtt [ rff w^e™*?^^ X 
Jo Jo Jo 



X 



t^x'^v'^ 



{e^^\P ■ xf + P^(e(^) ■ x){P ■ x)—^B^'\u) + {x'e^^^ - x^{e^^^ ■ x)) B^'\u) 



(6.47) 



Including the leading twist term, the total derivative term ( 6.44 ) and this last equation, 
equations for A^^^u), B^^^u), C(^)(m), A^^\u) and M"^ (u) can be obtained : 



u 



duB^'Hu)e'^'''' 



g{u) (6.48) 

f dug{u) I dte''^^'^ 
Jo 

{P-xf f du f dtt f dvv^e''"^^''B^'\u) 
Jo Jo Jo 
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1 1 



2(P-x)2 



(6.49) 

/'rfMC('=)(u)e'«^" = f du^^g{u) f dte f dvv'^e''"^'''' 
Jo Jo s Jo Jo 



duB^^\u) f dtt f dvv'^e'^'^^^ 
Jo Jo 



(6.50) 

rfM5(^)(M)e'«^" = [' dug{u){^^ + 1) f dtt^ f v^dve''"^'''' 
Jo Jo Jo Jo 

+ 4 / duB^^\u) f dtt f 
Jo Jo Jo 



{P-xf / duB^^\u) I dtt^ I dvv^e''^^^'' 

(6.51) 



1 







[' dug{u) (^2 + 1) f dt e f dv {z^{P ■ x) + 1) t2e»t?P- 
Jo Jo Jo 

+ 4 /" duB^^\u) [ dtt [ dvv'^e'"'^^'' 
Jo Jo Jo 



(6.52) 



Theoretically, all the chiral-even distributions B^^\u), C^^\u), A^^^ {u) can be computed 
once the leading (physical-) twist contribution A^'^\u) is known. But, practically, it quite 
long and hard if one has a conformal expansion for A^*^^ (u) (except if only the asymptotic 
part of the wave function is taken, see [|lll). The advantage of these equations is that they 
are independent of the model used to computed A^^'> (u) (another expansion can be used as 
the conformal one). They can be the basis of a numerical computation if one has discrete 
values for A^'^\u). 

In the appendix ^, similar equations for the chiral-odd distributions are given. 

6.3 Mass correction using the equation of motion 

Equation of motion is used to relate the two-point wave functions to the three points ones 
(these calculations come from [|Tl[]). 
From page 245, one has : 

idfju{x)'y^d{—x) = / dvu{x)gG n,y{vx)'y'^ d{—x) (6.53) 



The symbol 5^ is the derivative over the total translation define in p. 8 



25 



When the matrix element (0| . . . \p{P, A)) is taken (replacing x hj z), on gets (using the 
definitions [A.3| and |A.18| for (f)\\{u), gsiu), $(m) and "^{u)) : 



1 3eW-z 

2 " p-z Jo 



(6.54) 



= -ifpmlie^^'^ ■ z) j dv j Va[2<l>{a) + ^(a)] e-^P^(°"-"'^+''"s) 

The integration of the right hand side can be transformed : 

dv / PaF(a)e-'P"(""-"^+""«) = 2 dwe'^'^^"^-^^ / da^ / dau—F{a) (6.55) 
1 J Jo Jo Jo oig 

And a relation between the two-points and three-points distributions is obtained : 

I [011 (m) - g^iu)] [ dad I dau— [2<l>(«) + ^(a)] (6.56) 
In a similar way, the operator relation (in [^) 

u{x)'ypd{—x) = —i dv vu{x)x" gGai3{vx)^^d{—x) (6.57) 



gives 



-!('=)(«) + 2C^'\m) = / darf / cia^a^M - ad(l - m)) — [2$(a) + ^(a)] (6.58) 
8 Jo Jo «ff 

with 

C^'\u) = r r dtC^''\t) (6.59) 
Jo Jo 

Hence from the wave functions computed in section ^, the mass correction of the two- 
points functions can be computed. The advantage of this method is that it is quite simple, 
but it depends on the model for $ and and there is no general methods to find equations 
that relate the two-point to the three points functions. 

7 Conclusion 

In this work, I describe the different wave functions of meson, the method conformal 
expansion and its application to the non-leading twist wave functions. I also give methods 
of computing the mass correction, which helps to apply these results to heavier vector 
mesons as the p. 

The use of conformal expansion and exact operators relations coming from the equation 
of motion is still long and heavy, a more general method would be certainly useful, and 
there is no work of higher gluonic corrections, like four-points wave functions for example. 
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These non-local, non-perturbative objects are needed for the computation of exclusive 
decays, but they also describe the nature of the quark confinement in meson in a Lorentz- 
invariant and gauge invariant way. A better knowledge of the non-leading twist wave 
functions would probably help to understand how QCD describes the physics of hadron. 

Another question is the possible link to effective low energy theories; for example, the 
quark condensate (0|gg|0) is an order parameter of the breakdown of the chiral symmetry 
at low energy; can one found such thing with non-local wave functions ? 

The basis of the conformal expansion is the conformal invariance of QCD to one loop. 
But a analysis of the two-loop effects on this expansion, a better renormalization group 
dependence of these wave functions would be interesting. 

Other ways of determining these wave functions (by experiments or using another ex- 
pansion as the conformal one) is certainly desirable, and the "theoretical twist method" of 
determining the mass corrections that I give in this work will still be applicable. 
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A List of vector-meson wave functions up to twist 4 

In this appendix, 1 give a list of vector- meson wave functions (the chosen meson is the p). 



following the notation of P and fTl 



A.l Two-points distributions 

Chiral-odd p two-point distributions : 

{Q\u{z)a^,{z,-z\d{-z)\p-{P,\)) 



T 







2 ^1 



^ P ' ^ Jo 



{A.l] 



{0\uiz)[z,-z]di-z)\p-{P,X)) = -I (^/J - /p^^^^^ j (eW ■^)mj d«e^«^-/.(f)(«,p^), 

(A.2) 

Chiral-even two-points distributions : 

{Q\u{z)-fp[z,-z\d{-z)\p-{P,\)) = 
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Twist 


2 


3 


4 






0{l/Q) 














0± 


(v) (a) 
51 , 


/l3 



Table 1: "Physical twist" classification of two-points distributions 



P-z Jo 



1 e(^) ■ z 



2 '"{p-zy"^" 







(A.3) 



and 



\u{z)J^,-f5[z, -z]d{-z)\p- {P, \)) 
1 
2 



n \ jp Jp 



m. 



p^zp fdue'^P-^g'f\u,i,'). (A.4) 



with the notation 



^ = u- (l-u) = 2u-l. 



The vector and tensor decay constants fp and /J are defined as usually as 
(0|iZ(0)7^rf(0)|p-(P,A)) = /pm.eW, 



(A.5) 
(A.6) 



If every terms are needed, x'^ corrections has to be made. In that case, the notation 
is a little different : 
Chiral-odd two-points distributions : 

{0\u{x)ai^4^,-x]d{-x)\p~{P,X)) = 



X 



(P-X)2 P 

m 



1 / 2 2 



due'^^-'^C^^'Hu) 



with 



P-xJo 



P(°)(^) = /,(*)(^)-i0^(^,)-l/,3(^,) 



(A.7) 

(A.8) 
(A.9) 
(A.IO) 
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Chiral-even two-points distributions : 

{0\u{x)^^[x,-x]d{-x)\p-{P,X)) = 

,(A) ™ ^1 

fpmp 



^ P-x 



2 2 



X 



2 2 



1 e^-^) • a; 



— -a; 



2 ''{P-x)^ ''Jo 



with 



9± [u) 



u 



(A.ll) 



(A.12) 
(A.13) 
(A.14) 



A. 2 Three-points distributions 

Chiral-odd three-points distributions : 

{0\u{z)aap [z, vz\gG^,,{vz) [vz, -z\d{-z)\p~ (P, A)) 



T2 



/I 



e(^) ■ z 



2{p- Z) [^"^^^/^'^ ~ - Pc^Pudpt. + PpPi^Qai.Usp^p'^iv^Pz) 

+ \Pae^±l9pu - Ppe^±l9L - Pae^^lgp^+Ppe^^lga^\'^i{v,pz) 
+ [PMeia^/3z. -P^el^^az/ -P^^e^^^^^ +p^el^^^^]72(T;,p^) 

+ \Pae^_i^lp,,z^ - ppe'^llp^Za - Pae'ilpuZ^ + ppe^llpuZa]T4{v,pz) 

(A.15) 

where T is a twist 3 distribution ("physical twist"), the other are twist 4. There are other 
Lorentz-structures which are twist 5. 



{0\u{z)[z,vz]gGi^^{vz)[vz,-z]d{-z)\p {P, X)) = 



fp^l[(^±lpi^ - e±,,]S{v,pz) 



(A.16) 



{0\u{z)[z,vz]i-f5gGi^^{vz)[vz, -z]d{-z)\p (P, A)) 
■ fT 2r W 

These and S{v,pz) are twist 4 distributions. 



ifjmy^lp. - e^,]S(v,pz) (A.17) 
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Chiral-even three-points distributions : 

{0\u{z)i-fa[z,vz]gG^y{vz)[vz, -z]d{-z)\p' {P, A)) 
= fp^p 



■ z , 

+ [Pf.9au-Pu9ap]^i'^^Pz) 

- e(^) ■ z 

' {p- zf ^^^^" ~ P'^^f^P»'^('"^P^) 



(A.18) 



{0\u{z)ja-f5[z,vz]gGf,^{vz)[vz, -z]d{-z)\p (P, A)) 



Pa[Ppe^ll - p^e'^^l^jAiv^pz) 



e(^) ■ z 
p ■ z 
■ z 

' {p ■ zY 



[Pi.9L-Pu9ip]Hv,pz) 

[PpZy - PyZp]pa^{v,pz) 



where A{v,pz) and V{v,pz) are twist 3 distributions, the other twist 4. 
Every functions J-'{v,pz) have such form : 



(A.19) 



(A.20) 



and a is the set of three momentum fractions a = {a^, a^}. The integration measure 
is defined as 



/Va = / dad / dau / dag 6(1 — 
Jo io Jo 



ai . 



(A.21) 



B Models for p wave functions up to twist 4 



I reproduce the results of [|T^] here, which helps to understand section |]. Their computa- 
tions can be found in |§] and fll]], where a more complete list can be found. 



B.l Chiral-even distributions 

The model for the leading twist 2 distribution amplitude 0|j is 

3 



(f)\\{u) = 6uu 



l + 3a\^ + al^{5e-l] 



iB.l] 



with parameter values as specified in Tab. 0. The expressions for higher-twist distributions 
given below correspond to the simplest self-consistent approximation which satisfies the 
QCD equations of motion [H, 111 
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Three-particle distributions of twist 3 : 

V(«) = 540(3^1' 
A{a) = 360 Cattrfa^a^ 



540 (3^1 (ad - au)adaual, 



l + ^3^-(7«,-3)]. 



(B.2) 
(B.3) 



Two-particle distributions of twist 3 : 



(1) / \ 



1 + + <; + ^ Cs ( 1 



3 A 9 y 
16 ^ 16 ^ 



{^e - 1) 



+ 6 5+ {?>uu + -u In M + M In m) + 6 5_ (-u In -u — -u In m) , 



(B.4) 



4 
+ 



9 

m 



2 

15 
64 



+ 5C3) {^e - 1) 



C3 (3 UJ, 



(3 - 30^2 + 35^^) 



3 ~ 3 ~ 

+ -5+ (2 + Inw + Inu) + (2^ + Inw - Inu), 



(B.5) 



Three-particle distributions of twist 4 : 



$(a) 
$(a) 



3C3 + -C4 
3C3 + 3(4^ 



30(1 — ag)a 



5' 



30(a„ - ad)al 



120auadag, 



0. 



(B.6) 



Two-particle distributions of twist 4 : 



4 20 8 

5 + y^^+9^^ 



Qu{l — u) 



30^2(1 -m)^ 



10 20 

y C4 - y Cs 



■3 10 10/ 

2 + y^^ + y^^ 



(B.7) 



where the dimensionless couplings and (^4 are defined as local matrix elements 



(0|m5(G^^7c,75C?|p (P, A)) = /pm^Cs 



1 

3" 



+ ^/p"^JC4 



(B.8 



and have been estimated from QCD sum rules |T^. The terms in 5-1- and 5-1- specify 
quark-mass corrections in twist 3 distributions induced by the equations of motion. The 
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numerical values of these and other coefficients are listed in Tabs. H and |[|. Note that 
SU(3) breaking effects are neglected in twist 4 distributions and in gluonic parts of twist 3 
distributions. 

The numbers a\, a|, u^, 00'^ are computed in the appendix C of [Q, using usual two- 
points QCD sum rule. 5+ and 5_ have following definitions : 



/J m„ ± mg 
fp mp 



(B.9) 



B.2 Chiral-odd distributions 

The model for the leading twist 2 distribution amplitude 



IS 



Quu 



l + 3a^e + 4o(5e'-l) 



(B.IO) 



with parameter values as specified in Tab. 0. The expressions for higher-twist distributions 
given below correspond to the simplest self-consistent approximation which satisfies all 
QCD equations of motion [g, : 

• Three-particle distribution of twist 3 : 



r(a) = 540 C3UJJ {ad - au)adaual 



(B.ll) 



Two-particle distributions of twist 3 : 



h^;\u) 



Quu 



1 + 



4 



Cs^i (5^ 



-1-3 5+ {?>uu + ulnu + u\n.u) + ?>5^ {uXnu — u\n.u) 



(B.12) 



- 3e^ + ^a^e(3e^-l) + ^a^e^(5e^-3) + i^C3C^3^(3-30e^ + 35a 



hf{u) 



3 3 
+ -5+{l + i\Yiu/u) + -5^i{2 + \nu + \Yiu) 



(B.13) 



Three-particle distributions of twist 4 : 



0, 



S{a) 
S{a) 



-30Cj(ad - au)al 
30Cj(l — ag)a 

3ocr(i 



5' 



(B.14) 



Two-particle distributions of twist 4 : 

/isH = 6M(i-M)+5[cr+cr](i-3r 



u 



?,Ou^{l-u? 



2 4 T _ 8-7, 
5 + 3^^ 3^^ 



(B.15) 



^In the notations of Ref. 



ioj^o^ coi, Ci ^ G and ^ mm^I 
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V 








</> 


fxAMeV] 

J V L^vJ-c V J 


198 ± 7 


226 ± 28 


226 ± 28 


254 ± 3 


/^[MeV] 


160 ± 10 
152 ± 9 


185 ± 10 
175 ± 9 


185 ± 10 
175 ± 9 


215 ± 15 
204 ± 14 


aj 





0.19 ± 0.05 


-0.19 ± 0.05 





al 


0.18 ± 0.10 


0.06 ± 0.06 


0.06 ± 0.06 


0±0.1 




0.16 ±0.09 


0.05 ±0.05 


0.05 ±0.05 




ai 





0.20 ± 0.05 


-0.20 ± 0.05 







0.18 ±0.05 


-0.18 ±0.05 






0.20 ± 0.10 


0.04 ± 0.04 


0.04 ± 0.04 


0±0.1 




0.17 ±0.09 


0.03 ±0.03 


0.03 ±0.03 









0.24 
0.22 


0.24 
0.22 


0.46 
0.41 







-0.24 
-0.22 


0.24 
0.22 










0.16 
0.13 


0.16 
0.13 


0.33 
0.27 







-0.16 
-0.13 


0.16 
0.13 






Table 2: Masses and couplings of vector meson distribution amplitudes including SU(3) 
breaking. In cases that two values are given, the upper one corresponds to the scale 
^'^ — 1 GeV^ and the lower one to iJ? = m\ — ml = 4.8 GeV^, respectively, is taken at 
II — 1 GeV {ms — 150 MeV and the quark masses m„ and are put to zero). 





Cs 






T 
^3 


C4 


CI 


CI 




0.032 


-2.1 


3.8 


7.0 


0.15 


0.10 


-0.10 


V 


















0.023 


-1.8 


3.7 


7.5 


0.13 


0.07 


-0.07 



Table 3: Couplings for twist 3 and 4 distribution amplitudes for which we do not include 
SU(3) breaking. Renormalization scale as in previous table. 
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The constants (J and (J are defined as 

{0\ugG,.d\p~{P,X)) = 2/Jm^cr(eWp.-e(^)p^), 

{0\ugG,,ij,d\p-{P,\)) = ^fJmlCliel'^Pu-ei'^P,) (B.16) 

and have been estimated in from QCD sum rules: 

Cr~-Cj~0.10. (B.17) 

Other parameters are given in Tab. Like in the chiral-even case, SU(3) breaking cor- 
rections are neglected in twist 4 distributions. 

The number uj, a^, are computed in the appendix C of using usual two-points 
QCD sum rule. 5+ and 6- have the following definition : 

fpmu± rrid 



fi 

J o 



(B.18) 



C Formulas for Orthogonal Polynomials 

This appendix reproduces useful formulas for the Appell, Jacobi and Gegenbauer polyno- 
mials. 

Differentiation formula for Gegenbauer polynomials : 

d 



d^ 



Differentiation formulae for Jacobi polynomials : 



2 



(C.l) 



(C.2) 



'd^ 



(n + l) 



(n + 2)(2n + 3) 



(C.3) 



The equation ( |C.3| ) is obtained from ( |C.2| ) combined with ( p.4|) below. 
Recurrence formulae for Jacobi polynomials : 

(^ + l)(n + 3) (,,^^^^^,,,^^^^^ 



(n + 2)(2n + 3)' 



2n + 3 



Pi'^^iO 



+ 

n + 1 

2n + l 
n + l 

2n + l 



^ ^(1,0), 



2n + l 



^ -,(0,1), 



2n + l 



P^.'iO, 



(C.4) 
(C.5) 

(C.6) 
(C.7) 



Hn notations of Ref. MCl = (3/28)C3wJ. 
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Relations between Jacobi and Gegenbauer polynomials : 

(1 + o^rno + (1 - o^r^o = 2cy2(0, 

Orthogonality relations for Appell polynomials |T^ : 



Vaadaua Jk,i{ad,au)Jm,n{'^d,Oiu) = 6k+i,m+n 



(C.8) 
(C.9) 

(C.IO) 
(C.ll) 

(C.12) 

(k+l+l) 



where W, 



(k+l+l) 



k,m 



2*=+'+3(A; + Z + 3)(2A; + 2/ + 5)!!"*^'™ ' 

(C.13) 

d"^'^'^ Jk,i{ad, au) / da'^da'^ is a (k + l + l) x (A; + / + l) symmetric matrix. 



This result can be obtained from the following relations: 







(m + n < k + I) 



■2k+i+3(^k + l + 3){2k + 2l + 5)U 

while the integral is in general nonzero for m + n > k + I. 
Integral formulae for Appell polynomials : 



- {m + n = k + l), 



(C.14) 



d 



dad I da. 



du Jq Jo 1 — Oid — au \ oa, 



d d 

adT{ H auT: 1 ) - - au)^Jk,i{<^d, a«) 



dau 



uu k\l\{-lY 



k-l 



(1,1) 



2 {k + l + 2)\ \k + l + ?, 

d_ 

du 



(C.15) 



dad / dauadauil - - au)Jk,i{oid, a^) 
uu kW.{-l)'' ( k-l ^(1,1) 



2 (k + l + 3)\\k + l + 3^+'+^^^' n+ml^J 



(C.16) 



-f- I dad I dau 

du 



1 



d 



d 



1 — ttrf — \ attd oa^ 



adau{l - ad - auf JkA'^d-, 



Pl^:'UO + {k-l)Pl^[>{i) 



(1,1) 



fc!/!(-l)fe r (A; + / + 2)(fc + / + 4) 

'J {k + 1 + 3)1 [ k + l + 3 ' k+l+2\^l TV- - V- fc+i+lV 

dad I dau 



(C.17) 



1 



I - ad - au V "^Od 9a„ 



9 d 

+ 



adau{l - ad - au)^Jk,i{'^d, ««) 



kW.{-l)'' 
{k + l + 3)14 



-k + l 



2k + 2l + 5\ ^■■+'+^ 



k + l + 2 
2k + 21 + 



_ ( p(0,0) _ p{0,0) /px 



(C.18) 
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The results (|C.15|) — (|C.17|) can be obtained by differentiating and/or integrating the Appell 
polynomials Jk,i{c(d, c^u) term by term. To obtain ( |C18| ), it is convenient to calculate its 
derivative first, which can be done similarly to ( |C.15|) , (|C.16|) and ( |C.17|) , and then integrate 
the result with the condition that it vanishes at m = 0. 



D Operators relations for different theoretical twist 



This appendix contains some relations taken from |T^, and which help to isolate 
the different twist parts of two-points non-local operators. 

D.l Chiral-even operator 

Consider the non-local operators ilj{x)'ja'^{—x). The local expansion around x = can be 
written in the following way : 

oo ^ 

^{x)-fa^{-x) = • ..Xf^^—tfj D^, ... Df,„ la^ (D.l) 

n=0 

In order to get the leading twist contribution, the symmetrization has to be done over 
all indices and subtract the traces : 

CO ^ 

[?(a;)7aV^(-a;)] 1 t =y2x^^...x^—i) ( B^i • • • 7a ) - traces (D.2) 

n\ \ / sym 

71=0 

The symmetrization has an integrated solution : 

V'(x)7a^(-X) = 
J sym 

V— -\ x^^ x^" -fl^ n ^ ^ ^ 1 

= H ^(0) 1 ^ D,,... 7. + ^ DaD,, . . . D,^_, 7m„ [^(0) • 

n=0 ' ^ 

(D.3) 

' d - 

dv — — ■iIj{vx)xi/j{—vx) (D.4) 







dxa 

with X = '~f^x'^. 

In order to have really the leading twist part of ilj{x)'ja'4^{—x), the leading twist part of 
■il!{x)xilj{—x) is needed. If one does a local expansion of this latter non-local operator, one 
will already have symetrized local operators. But traces has alos to be removed; this can 
be expressed as a differential equation : 

^ = (D.5) 

which has the formal solution : 

^{x)xip{—x)~\^ ^ = tjj{x)xijj{—x) 

(D.6) 
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Equations for the twist 3 part of ■ip{x)'ya4'{—x) are the following 



Qs I du dvip{ux) uG^,y{vx)x'^x^5 



twist 3 JQ 

— ivG^u{vx)x'^x ip{—ux) 



^(a;)7M75V^(-a;) 



twist 3 



?/'(mx)7/375V'(-mx) 



du I dvip{ux) uG^u{vx)x^x 



— wG^^{vx)x'^x'-y^ tp{—ux) 



+ ie 



uaf3 



udu x,,dn 



ip{ux)'~^pip{—ux) 



(D.7) 



where G^y is the gluon field strength, G^y = {l/2)e^uapG°'^ , and da is the derivative over 
the total translation : 



da 



tp{ux)'Ji3lp{—Ux) 



_d_ 

dye 



ip{ux + y)ji3ilj{—ux + y) 



(D.8) 



These equations come from and can be obtained with the equation of motion, 
neglecting quark masses. 



Another useful formula is the equation (5.13) in ||T9| which can be written like this : 

(D.9) 



d — ^ — . 

'?/^(a;)x-?/'(— x) = —d'^ip{x)xilj{—x) + operators with gluons 



dxadx" 

The expansion of p.6| at the order x^, gives 

1 



ilj{x)x^jj{—x) 



1 



■ip{vx)xijj{—vx) + . . . 



[ijj{x)xijj{—x)~\^ ^ — -x^ I dvv'^d^ip{vx)x^jj{—vx) 



JO 



+ operators with gluons + twist 3 + . . . (D.IO) 
This last equation gives a twist 4 contribution with a total derivative (defined in p.8|) . 

D.2 Chiral-odd operator 

Consider the non-local operators '?/^(x)crQ^'?/'(— x). To get the most symmetrical part in the 
local expansion, one can compare it with the solution in p.4| and get directly : 



ij{x)aal3lp{-x) = 
J sym 

1 d - 

dv — — ip{vx)aaux''ip{—vx) — {a ^ P) 







dx 



(D.ii; 



/3 



37 



In 0, the equations for the leading twist part of ip{x)aaux"''p{—x) can be founds : 

d 



dXa 

dxpdxP 



[ilj{x)aaux''ip{-x)] 



1. t. 











(D.12) 
(D.13) 



Like in equation p.lO| , the contribution of total derivative can be added : 



4j{x)aaux''ip{-x) 



1 



c)aaux''ip{-x)\ - 
+ operators with gluons + . . . 



dtt^d'^ {^{tx)aau4^{-tx)) 

(D.14) 



E Mass corrections to two-points distributions 

In this appendix, the mass corrections for the distributions are given, using the method 
described in section |0. 



E.l Chiral even distributions 

For the different wave functions 

{0\u{x)jp[x,-x]d{~x)\p-iP,\)) = 



P„ 



P-x 



1 1 
X m. 



4 



+ eW - P 



g(A) . ^ 



X 



2 2 

X mt 



1 pW . T 



one has 



2 '^{P-x] 



duB^'\u)e'^^'' 







1 1 



2 {P-x) 



duC^'\u)e 



iiPx 



f duA{e){u) [ dte'*«-P^ 
Jo 

-{P-xf f du f dtt [ dw^e^'^^^'^B^^h 
Jo Jo Jo 



du^-^A^'Hu) ! dte ! dw^e'^'^P^ 



(E.l) 



u) 



(E.2) 



duB'-'\u)e'^^'' 



^ JO JO 

- [ duB^^\u) [ dtt [ dvv'^e''"^^'' 
Jo Jo Jo 



C duA^'^ {u) (e' + 1) f dt t^ f v^dve''''^^'' 
Jo Jo 



(E.3) 
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4 [ duB^^\u) [ dtt [ dvv 
Jo Jo Jo 

{P-xf ! duB^^\u) f dtt^ [ dvv^e^ 
Jo Jo Jo 



itv^Px 

(E.4) 



[ duA^''\u)e'^^'' ^ 
Jo 

f duA^'\u){^'' + 1) f dtt" f dvv^ {i^{P ■x) + l) i2e*''*«-P- 
Jo Jo Jo 

+ 4 [ duB^^\u) f dtt j dvv^e'"^^^'' 
Jo Jo Jo 



E.2 Chiral odd distributions 

For the different wave functions 

{Q\u{x)a^,^[x, -x\d{-x)\p~ {P,X)) = 



~l~ {,Pfj,Xv PvXfi) 



e(^) • X 

2 ,-1 



-ra. 



one has 



Jo 



due'^^''A^''\u) = 

/ duE{u) / dt I dvv 
Ju Jo Jo 

lit 



S^ivt^Px 



+ 2- ^ 



tv^{Pxf v^{Pxf v^{Px)^ 



+2 



v'^{Pxy v{Px) v{Px) 



v{Px) 



Jo 



\{P-x) [ duE{u) f dt f dvv 
2 Jo Jo Jo 



2^itviPx 



(E.5) 



(E.6) 



(E.7) 



^ -4.-!$^-2(iC)H + 2t 



tv'^{Px) v{Px) 



(E.8) 
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v{Px)^ 2 ^^^^ 2 



tv^{Px) 



1 



+ 



-2(Px) 



1 



+ 



v{Px) 



(E.9) 



where 



E(u) = -t-^A^°\u) 



(E.IO) 



There is no equation for I3^^\u) because it is aheady a O(m^) contribution. 
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